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ABSTRACT 

/ 

The development of vehicles powered by direct-current 
sources together with the development of the silicon-control- 
led rectifier has led to the use of squirrel-cage induction 
motors operating on non-sinusoidal , variable-frequency volt- 
age supplies. A digital computer simulation of the tran- 
sient and steady-state pe r^o r^hW^o f a three-phase motor 
is derived and its use in predicting motor operation is 
illustrated. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 
(All rotor parameters referred to the stator) 

= Instantaneous value of voltage, volts 

= Instantaneous value of current, amoeres 

= Stator phase self-inductance, henrv (assumed 
identical for all three phases) 

= Rotor phase self-inductance, henry (assumed 
identical for all three phases) 

= Stator-rotor phase mutual inductance maximum, henry 

= Rotor position angle, radians 
= Rotor phase resistance ohms 
= Stator phase resistance, ohms 

= Electromechanical torque developed, newton-meters 
= Mechanical torque, newton-meters 

= Rotor angular velocity, radians/second 
= Inertia of system referred to motor, kilogram meters 
= Friction and windage coefficient 
= Magnetic flux linkages , weber 
= d/dt 

= Stator phases 
= Rotor phases 
= Time, seconds 
= Time increment, seconds 
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1 . 



INTRODUCTION 



Engineers have long sought a method of combining the 
simplicity, compactness, low cost, ruggedness and low main- 
tenance characteristics of the squirrel-cage rotor type of 
induction motor with the speed control characteristics of 
direct-current motors, involving comutators and brushes, or 
of wound- rotor induction motors also using brushes. 

In recent years two developments have given importance 
and feasibility to this project. 

1) The present and projected development of space ve- 
hicles, deep-submergence-oceanographic vessels, and nuclear 
reactor coolant pumps has demonstrated and intensified the 
requirement for reliable, variable-speed drives which reauire 
no maintenance, not even brush replacement, during their ex- 
pected lifetime, and which are very reliable. Preferably 
these drives should be able to operate from a direct-current 
source such as a fuel cell or battery pack. 

2) The development of the silicon-controlled rectifier 
and its use in bridge-inverters controlled by solid-state 
logic devices provides a reliable, practical method of ob- 
taining a variable-frequency three-phase power supply from 

a direct-current source. These inverters commonly contain 
no moving parts, have very short switching times and appear 
to be very reliable. 

Many descriptions of solid-state inverters have been 

r 9 1 

published in recent years and it will be assumed that 

voltage wave shapes of the following type are readily avail- 
able . 
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Power Conditioner Voltage Output 
Figure 1-1 



2. ASSUMPTIONS 

1) The wave shapes of all the three phase voltages are 
cyclic and have the same shape. 

2) Iron losses are assumed to be small and are not 
accounted for. Saturation of the magnetic path is neglected. 

3) Deep-bar effects are not considered. 

4) Since the source voltage is applied to one phase in 
series with two phases in parallel (Y connection) , and the 
windings are assumed to be symmetrical, the phase voltage 
division is one-third and two-thirds of the source voltage. 

5) The rotor is cylindrical as is the inside of the 
stator. The air gap is therefore uniform. 

6) The rotor is of the symmetrical, three-phase, coil- 
wound type. This is done to facilitate analysis although 
the basic phenomena occurring in a squirrel-cage rotor are 
the s ame . 



10 



3. 



PRINCIPLES OF OPERATION 




Single-Phase Motor Model 
Figure 3-1 

If at time zero a positive voltage v is applied to 
stator phase A as shown, current i A will flow as indicated 
creating flux which will link the turns of phase a. As 
increases an induced current will be generated in rotor 
phase a which will oppose the time rate of change of i^. 
Fleming's right-hand rule may be used to establish the po- 
larity of this induced current as negative, i.e., in at the 
bottom and out at the top. By Lenz's law a force or torque 
is created which acts to turn the rotor in a counterclock- 
wise direction. Consider the sequence of events: v A creates 

an exponentially rising current i A , i A produces flux linkage 
¥ , the time rate of change of which generates voltage v , 

and in turn v creates current i . 

a 

It is obvious that if periodic wave forms are consid- 
ered i a must lag i A . Also, since v^ is created by the time 
rate of change of i A , when i A begins to decrease in value 
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the polarity of induced voltage v a must reverse. Therefore, 

at this point it may be assumed that the wave form of i a 

must lag the wave form of i A « 

If now the rotor is considered free to rotate and it 

is assumed that counterclockwise torque is desired, the 

following relationships must be considered. Let 0, rotor 

angle, be zero when coil axes A and a are colinear and let 

counterclockwise be considered positive. Then the mutual 

inductance will equal L^ x cos 0 by definition. Torque, T, 

(3. 

will then be equal to i A * (-i a ) * tgg cos 0] which equals 



"V (_i a ) * [L lx sin 01 ‘ 

For positive torque with 0 greater than zero and less 
than II radians, the product i i must be negative. For 0 

A cl 

between II and 211 radians i. i must be positive. 

A cl 

Another effect which must be considered is the rotor 
phase belt "carrying" its current around with it. In the 
numerical methods solution to follow it will be assumed that 
electrical transients are fast enough compared to mechanical 
transients so that 0 may be considered constant for the very 
short At periods involved. The foregoing does not mean that 
"speed voltages" or counter-electromotive force is neglected, 
but only that p0 and mutual inductance are held constant for 
the At periods . 



12 



4. 



ANALYSIS OF THREE-PHASE MOTOR OPERATION 
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Figure 4-1 



If each coil is considered to have an axis along which 
a composite flux is created the total flux linkages along 
each axis may be computed as follows: 

\ m Vll + l B L yx cos (2n/3) + i C L yx cos <- 2n / 3 > + Vyx 

cos (0) + i, L cos (0+2JI/3) + i L cos (0-2JI/3) 
b yx ' c yx 

f B = i B L ll + i C L yx COS (21T/3) + i A L yx COS (_2II/3) + i a L yx 

cos(0-2II/3) + i.L cos (0) + i L cos(0+2II/3) 
d yx c yx 

31 C " i C L U + V-yx cos (2n/3) + i B L yx cos t-2n/3) + i a L yx 

cos(0+2II/3) + i,L cos (0-211/3) + i L cos (0) 

' b yx c yx 

‘'a = ^xx + i b L yx cos (2n/3) + i c L yx cos ( ' 2n/3 + i A L yx 
cos (-0) + i^L cos (2JI/3-0) + ipL cos ( — 211/ 3—0) 

f b = i b L xx + i a L yx COS ( ~ 2II/3) + L c L yx COS (2II/3) + i A L yx 

cos ( — 211/3—0) + i L cos (-0) + i n L cos (2JI/3-0) 

d yx ^ yx 
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1 'c * i o L xx + SV cos (2n/3) + i b L yx cos <- 2 ”/ 3) + iftV 

cos ( 2 n/3— 0 ) + i D L cos ( — 2IT/3— 0 ) + i^L cos (-0) 

B yx C yx 

The voltage equations for the physical machine are de- 
fined to be: 

v = iR + pT 

i . e. : 

V A = £ A R + X 11 pi A + L yx COS (2II/3) pi B 

+ L yx cos ( — 211/3) pi c + L^ x cos (0) pi a - i a L yx sin 

(0) p0 + L^ x cos (0+2II/3) pi b - i^L^ x sin (0+2II/3) p0 

+ L^ x cos ( 0-2II/3 ) pi c - i c L yx sin (0—211/3) p0 

v B = igR + L n pi B + L yx cos ( 2II/3) pi c + L yx cos ( — 211/3) 

pi A + L yx cos (0-2II/3) pi a - i a L yx sin (0-2II/3) p0 

+ L yx c °s (0) pi b - i b L yx sin (0) p0 + L yx cos 

(0+2II/3) pi - i L sin (0+2II/3) p0 

O U y 

V C = i c R + L n pi c + L yx cos (2II/3) pi A + L yx cos 

(-2K/3) pig + L yx cos (0+2II/3) pi & - i & L yx sin 

(0+2n/3) p0 + L cos ( 0-2II/3 ) pi, - i, L sin 
yx d D y x 

(0-2II/3) p0 + L cos (0) pi - i L sin (0) p0 
yx c c yx 

v a = i a r + L xx pi a + L yx cos (2II/3) pi b + L yx COS ( “ 2n / 3 > 
pi c + L y x C ° S (_9) pi A + i A L yx Sin (-0) pQ 
+ L yx C ° S (2II / 3_0 ) Pi- B + 1 B L yx Sin (2n / 3-0 ) P0 
+ L yx cos ( -2n / 3-0 ) pi c + i c L yx Sin (-2II/3-0) p0 
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v^ 



b = V + L XX ^b + L yx COS (- 2n / 3 ) Pi a + L yx COS (2n/3) 
pi + L cos (-2II/3-0) pi + i L sin (-211/3-0) p0 



yx 



A yx 



+ L cos (-0) pi~ + i n L sin (-0) p0 + L cos 
yx r B B yx ^ yx 

(2n/3-0) pi c + i c L yx sin (2n/3-0) p0 

v c " V + L xx pi c + V cos (2n/3) pi a + L yx cos ( ‘ 2n/3) 
pi b + L yx cos (211/3-0) pi A + i A L yx sin (211/3-0) p0 

+ L yx cos (-2II/3-0) pig + igL yx sin (-2II/3-0) p0 

+ L yx cos (-0) pi c + i c L yx sin (-0) p0 

If, for simplicity of notation, terms are redefined as 

follows; D = L cos (2II/3) 

yx 

E = L cos ( — 211/3) 
yx ' 

F = L cos (2II/3-0) 
yx 

G = L cos (2II/3+0) 

yx 

H = L yx cos (-2H/3-0) 

I = L cos (-211/3+0) 

yx 

J = L cos (0) 

yx 

K = L cos (-0) 

yx 

M = L sin (0) p0 

yx 

N = L sin (-0) p0 

yx 

P = L yx sin (2II/3+0) p0 

Q = L sin (2TI/3— 0) p0 
yx 

S = L sin (-2II/3+0) p0 

yx 

T = L sin ( — 2II/3—©) p0 
yx 



15 



The equations in matrix form are: 
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0 




x a 


v b 
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0 




A b 


v c 
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0 


r 




_ 1 c_ 



L 11 D E J G I 




P X A 


E L i;l D I J G 




pi B 


D E L i;l G I J 




pi c 


K F H L xx D E 




P £ a 


H K F E L xx D 




P A b 


_F H K D E L xx 




_ pi c 



It may be well to note at this point that no assump- 
tions or specifications have been made as to voltage wave- 
form or frequency. It has been assumed that the motor wind- 
ings are symmetrical, an assumption which simplified the 
notation and in addition allowed the use of the original 
voltage distribution assumption for a three-wire system. 

The modification required to consider a four-wire unbalanced 
system is obvious and will not be pursued. 
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The system of nonlinear differential equations obtained 
may be solved only if 0 and p0 may be considered constant. 
Since the transient solution is desired this can only be -jus- 
tified for very short time periods. 

The application of numerical methods and the use of a 
computer may be expected to produce an inexact but useful 
solution to the problem. 

Let the previously derived system of matrix equations 
be expressed by: 

[v] = [H] [i] + [K] [pi] 

then : 

[pi] = [K] _1 [v] - [K ] _1 [H ] [i] . 

Torque generated in the constant air-gap machine is 
equal to the summation of the product of stator phase cur- 
rents, rotor phase currents, and the space rate of change in 
their respective mutual inductances. 

Therefore : 

T = -i fl i a L yx sin 0 - i A i b L yx sin (0+2II/3) - i A i c L yx 

sin (0-2II/3) - i D i,L sin 0 - i D i L sin (0+2IT/3) 

d jo yx Jj c yx 

- VaV Sin <®- 2n / 3 > - Vc L yx sin 0 ‘ VaV 

sin (0+2JI/3) - i i, L sin (0-2II/3) 

w Jo y x 

The mechanical equation for torque: 

T - T m = J(p(p0) ) + Zp0 

may be solved for p(p0): 

p(p0) = (T-T m -Zp0)/J . 
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If load torque is assumed to be known, predictor- 
corrector methods may be used as follows: 

PREDICTOR 

(pi) n = f(0,pO,i) n 
p(p0) n = f ( i ,0) n 
W = i n + <*t) (pi) n 



CORRECTOR 



Check 



(pi) n = f (0/p0/i) n+1 
I n+1 = I n + dt (pI+PI)/2 



Error 




n+1 



If the absolute value of error is found to be unaccept- 
able, set I n+ j_ = an d recompute I_ n+ ^ . If the error val- 

ues are acceptable, torques, rotor angles and angular veloc- 
ity may be determined preparatory to the next step. Adjust- 
ing At after each cycle of computations will permit saving 
of computer time while maintaining the desired degree of 
accuracy. 

Note that all motor parameters used in determining ex- 
pected transients or steady-state response may be determined 

. f 5 1 

by relatively simple laboratory tests of a given machine. 

2 

Stator and rotor copper (I R) loss calculations may be in- 
corporated into the program. 

It is theoretically possible to account for hysteresis 
and eddy-current losses by detailed analysis of flux changes 
and the magnetic path. However, the parameters required for 
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these calculations are not normally available to the engineer 
for an "off the shelf" motor. Since the induction motor 
circuits themselves are heavily inductive it may ontimis- 
tically be assumed that the circuits act as a low-pass filter 
and typical loss figures for the same class of motor opera- 
ting on a three-phase balanced sinusoidal voltage suooly may 
be used. 

5. APPLICATION 

For purposes of illustration a 20-horsepower, 440-volt, 
three-phase, 60-Hz., 3470-rpm, two-oole, squirrel-cage in- 
duction motor is selected. The stator windings are WYE con- 
nected and the phase model parameters, as determined by 
blocked-rotor and no-load tests ^ are as follows: 

R = 0.30 ohms 
r = 0.36 ohms 

= 0.00167 henrys (leakage inductances) 
model = 0.0723 henrys (magnetizing inductances) 

The model magnetizing inductance being three-halves the 
actual value, the inductance values used are: 

L 11 = L i + L yx (2 / 3) = 0-04987 

L = L, , = 0.04987 
xx 11 

L = L model (2/3) = 0.0482 
yx yx 

The mechanical load is arbitrarily assumed to be a pro- 
peller with load torque equal to 0.000315 times the square 
of the angular velocity in radians per second. This coeffi- 
cient is selected to require 19.5 horsepower delivered to 
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the load at rated speed and the friction coefficient (0.0282) 
accounts for the remaining one-half horsepower at rated con- 
ditions . 

The program following, Appendix I, was successfully run 
on an IBM 360 computer. Limiting the storage requirements 
to that required for 200 steps at a time allowed double- 
precision computation using less than 100,000 bytes of core. 

IBM supplied the subroutine for matrix inversion (GAUSS3) . 
Provision has been made to avoid multiplication or division 
by near-zero or zero values. 

SINUSOIDAL VOLTAGE SUPPLY 

In order to obtain a standard of comparison, a three- 
phase, 60-Hz, sinusiodal voltage supply was assumed for simu- 
lation purposes. Figures 6-1 through 6-6 display the com- 
puted current wave-forms and generated torque. The areas of 
special interest are considered to be motor starting tran- 
sients and steady-state performance. 

NON -SINUSOIDAL VOLTAGE SUPPLY 

A 60-Hz, pulsed voltage supply is assumed as an output 
from the power conditioner which, in accordance with assump- 
tion number four, is resolved in the ratio one- third and 
two-thirds across the stator phases. Computed phase voltage 
wave-forms, current wave-forms and generated torque values 
are plotted against time for starting and steady-state con- 
ditions in Figures 6-7 through 6-10. Voltage magnitudes are 
computed to have the same rms value as the rated sinusoidal 
voltage . 
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VARIABLE -FREQUENCY VOLTAGE SUPPLY 

Assumptions are the same as above except that tachom- 
eter feedback is assumed to control the frequency and thus 
the pulse widths of the stator voltage supply. If the cir- 
cuit model Thevenin equivalent is solved as a function of 
phase-voltage frequency and rotor angular velocity for maxi- 
mum power transfer to the load, and thus maximum generated 
torque, and if stator-phase resistance is neglected, the op- 
timum voltage frequency is found to be equal to rotor angu- 
lar velocity plus a constant. The results of this simulation 
are recorded in Figures 6-11 through 6-13. 

While very large starting torques are obtained, the 
currents generated are excessive and any motor not especially 
designed would probably be destroyed by a short period of 
operation involving a series of mechanical transients. The 
current requirements and heat dissipation requirements for 
the power conditioner would also far exceed those required 
for steady-state operation. 

Consider the steady-state relations that magnetic flux 
linkage per pole is proportional to the ratio of rms voltaae 
to frequency and that torque generated is proportional to 
the space rate of change of magnetic flux linkages. ^ Con- 
sider also that the variable-frequency pulse width described 
above is one-sixth of the stator- voltage cycle time. The 
rms voltage for a given frequency may be controlled by vary- 
ing the pulse width up to a maximum of one-sixth of the cycle 
time. If the pulse width is made proportional to f requencv , 
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the conditions quoted above are satisfied; therefore, maximum 
rated torque will be available during mechanical transients, 
and rms currents will be reduced. The results of this simula- 
tion, including the described voltage wave-form, are recorded 
as Figures 6-14 through 6-17. 

6. CONCLUSIONS 

Results displayed in Figures 6-1 through 6-13 compare 

f 8 1 

favorably with those obtained by other methods . Figure 

6-18 is a comparison of p© vs. time for a given motor opera- 
ting on different voltage wave-forms. 

From the data obtained the following conclusions may be 
drawn : 

The pulsed wave-form output obtainable from a power 
conditioner is useable for squirrel-cage motor 
operation. 

Variable speed control may be obtained by variation 
of power conditioner frequency. 

If the load requirements are such that few large 
speed changes are expected, such as normal opera- 
tion of a nuclear reactor coolant pump , the added 
complexity of the power conditioner required to con- 
trol pulse-width will probably not be justified. 

If large and frequent speed changes are expected, 
as would probably be the case in propelling an ocean- 
ographic vessel, the added complexity involved in 
pulse-width control is justified. 
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The computer program used to simulate the last wave- form 
described is presented as Appendix I. Computations using 60 
time increments per stator voltage cycle required approx- 
imately 11 minutes of IBM 360 Computer time to simulate two 
and one-half seconds of motor transient time. 

Areas of further investigation should include considera- 
tion of the effects of magnetic saturation and its effect on 
current and generated torque. Inclusion of power conditioner 
operation, copper losses, and vibration calculations in the 
simulation would improve the program's validity and thus its 
usefulness . 
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Stator Phase Currents vs. Time. Sinusoidal 
Voltage Supply. Motor Undergoing Starting 
Transients. 



Figure 6-1 



24 




Rotor Phase Currents vs. Time. Sinusoidal Volt- 
age Supply. Motor Undergoing Starting Tran- 
sients. 



Figure 6-2 
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Torque Generated vs. Time. Sinusoidal Voltage 
Supply. Motor Undergoing Starting Transients. 

Figure 6-3 
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Stator Phase Currents vs. Time. Sinusoidal 
Voltage Supply. Motor in Steady-State Opera- 
tion. 



Figure 6-4 
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Rotor Phase Currents vs. Time. Sinusoidal 
Voltage Supply. Motor in Steady-State Opera- 
tion. 



Figure 6-5 
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Torque Generated vs. Time. Sinusoidal Voltage 
Supply. Motor in Steady-State Operation. 

Figure 6-6 
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Stator Phase Voltages vs. Time. Constant - 
Frequency, Non-Sinueoidal Voltage Supply. 
Motor Undergoing Starting Transients. 

Figure 6-7 
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Stator Phase Currents vs. Time, Constant- 
Frequency, Non- Sinusoidal Voltage Supply. 
Motor Undergoing Starting Transients. 

Figure 6-8 
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Rotor Phase Currents vs. Time. Constant- 
Frequency, Non-Sinusoidal Voltage Supply. 
Motor Undergoing Starting Transients. 

Figure 6-9 
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Generated Torque vs. Time. Constant- Fre- 
quency, Non-Sinusoidal Voltage Supply. Motor 
Undergoing Starting Transients. 

Figure 6-10 
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Stator Phase Currents vs. Time. Non-Sinus oidal 
Voltage Supply With Frequency 90 Radians Per 
Second Greater Than Rotor Angular Velocity. 

Figure 6-11 
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Rotor Phase Currents V3. Time. Non- Sinusoidal 
Voltage Supply With Fi squency 90 Radians Per 
Second Greater Than Rotor Angular Velocity. 
Motor Undergoing Starting Transients. 

Figure 6-12 
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Generated Torque vs. Time. Non-SInusoldal Volt- 
age Supply With Frequency 90 Radians Per Second 
Greater Than Rotor Angular Velocity. Motor 
Undergoing Starting Transients. 

Figure 6-13 
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Stator Phase Voltages vs. Time. Variable Fre- 
quency, Variable Pulse-Width Non-Slnusoldal 
Voltage Supply. Motor Undergoing starting Tran- 
sients. 



Figure 6-14 
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Stator Phase Currents vs. Time. Variable Fre- 
quency, Variable Pulse-Width Non-Sinusoidal 
Voltage Supply. Motor Undergoing Starting Tran- 
sients. 



Figure 6-15 
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Rotor Phase Currents vs. Time. Variable Fre- 
quency, Variable Pulse-Width Non-Sinusoidal 
Voltage Supply. Motor Undergoing Starting Tran 
sients. 



Figure 6-16 
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Generated Torque vs. Time. Variable Frequency, 
Variable Pulse-Width Non- Sinusoidal Voltage 
Supply. Motor Undergoing Starting Transients. 

Figure 6-17 
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Rotor Angular Velocity vs. Time. Comparison of 
Motor Response to Different Voltage Inputs. 

Figure 6-18 
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APPENDIX I 



FORTRAN IV COMPUTER PROGRAM FOR SIMULATION OF THREE PHASE 
INDUCTION MOTOR PERFORMANCE ./HEN OPERATING ON THREE-PHASE^ VARIABLE 
PULSE-. /IDTH, VARIABLE FREQUENCY, SQUARE NAVE VOLTAGE SUPPLY: 
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VMAX = 540 



VSM = VM<s x/ A. 
VLR=VMAX-VSM 
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